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We study the low-energy limit of a quarter-filled one-dimensional Mott insulator. We analytically 
determine the local density of states in the presence of a strong impurity potential, which is modeled 
by a boundary. To this end we calculate the Green function using field theoretical methods. The 
Fourier transform of the local density of states shows signatures of a pinning of the spin-density 
wave at the impurity as well as several dispersing features at frequencies above the charge gap. 
These features can be interpreted as propagating spin and charge degrees of freedom. Their rela- 
tive strength can be attributed to the "quasi-fermionic" behavior of charge excitations with equal 
momenta. Furthermore, we discuss the effect of bound states localized at the impurity. Finally, 
we give an overview of the local density of states in various one-dimensional systems and discuss 
implications for scanning tunneling microscopy experiments. 

PACS numbers: 68.37.Ef, 71.10.Pm, 72.80.Sk 



I. INTRODUCTION 



Over the past decades scanning tunneling microscopy 
(STM) and spectroscopy techniques have been estab- 
lished as an important experimental tool to study 
strongly correlated electron systems with a spatial res- 
olution down to the atomic scale. In these experiments 
the tunneling current between the sample and the STM 
tip is measured as a function of its position and the ap- 
plied voltage. The tunneling current is directly related^ 
to the local density of states (LDOS) in the sample; thus 
STM experiments provide a tool to investigate the spa- 
tial dependence of the LDOS, which may arise for exam- 
ple due to impurities. These spatial modulations can be 
analyzed in terms of the Fourier transform of the tun- 
neling conductance which is directly proportional to the 
Fourier transform of the LDOS. Using this method of 
analyzing STM data one can infer informations about 
the bulk state of matter like the properties of its quasi- 
particle excitations. For example, this has been used to 
study high-temperature superconductors^ as well as car- 
bon nanotubes^. 

Theoretical studies of the LDOS and STM in the pres- 
ence of boundaries have been carried out in particular for 
one-dimensional systems including Luttinger liquids^"^, 
open Hubbard chains-'^ , and charge-density wave (CDW) 
statesiSiiii. In these systems the coupling to impurities 
is relevantiir— , which leads, at sufficiently low energies, 
to an effective cutting of the system into disconnected 
pieces. In this sense the effect of an impurity can be 
modeled by a boundary condition. Previous studies of 
the Fourier transform of the LDOS in both gapless^^- 
and gappediSiiii systems revealed a pinning of the CDW 
at the impurity as well as individually propagating spin 
and charge degrees of freedom. 

Here we study a quarter-filled one-dimensional Mott 
insulator in the presence of a strong impurity poten- 



tial. It has been establishedi^ii^ that for sufficiently 
strong interactions double Umklapp processes are rele- 
vant and generate a gap in the charge sector. As a mi- 
croscopic realization one may think of a quarter-filled 
Hubbard model extended by including nearest-neighbor 
interactions, which possessesi^ an insulating phase for 
sufficiently strong repulsions. Experimentally these sys- 
tems are relevant for the description of organic quasi-one- 
dimensional conductors, i.e. the Bechgaard and Fabre 
salts.^''. 

In this article we study the low-energy behavior of 
the LDOS by employing the field theoretical description 
of the Mott insulator. The spectrum consists of mass- 
less spin excitations (spinous) carrying spin ±1/2 but 
no charge and massive charge excitations (solitons and 
antisolitons) carrying charge ±e/2 but no spin. In this 
framework the impurity is modeled by a boundary for the 
collective excitations. In particular, we will focus on sig- 
natures of dispersing quasiparticles in the spatial Fourier 
transform of the LDOS, as they have been clearly ob- 
served in similar studies of Luttinger liquids^"— and half- 
filled Mott insulators as well as CDW statesi^iii. The 
spectral function of the translationally invariant, quarter- 
filled Mott insulator has been previously derivedii by 
Essler and Tsvelik. From the quantum numbers stated 
above it is clear that an electron will fractionalize into 
at least one spinon and two antisolitons. Thus the spec- 
tral function exhibits a featureless scattering continuum 
of at least three excitations, which is also consistent with 
experimental results^^ from angle-resolved photoelectron 
spectroscopy on Bechgaard salts. 

This article is organized as follows. In Sec. IH] we 
will briefly discuss the field theoretical description of the 
quarter-filled Mott insulator. In Sec. IIIII we present our 
results for the single-particle Green function. This will 
then be used to derive an analytic expression for the 
LDOS, which we will discuss in detail in Sec. IIVI Our 
main result is presented in Eq. (jl5p as well as Figs. [T] 



and [21 which show the Fourier transform of the LDOS in 
the presence of an impurity. Surprisingly, we find two dis- 
persing features which follow the dispersion relations (fT9|) 
and (pni) respectively. These features can be interpreted 
as arising from an antisoliton or a spinon-antisoliton pair 
propagating between the position of the STM tip and 
the impurity. In addition, we observe a non-dispersing 
singularity at momentum Q = 2kp which is indicative of 
a pinning of a spin-density wave (SDW) at the bound- 
ary. A detailed analysis further reveals dispersing fea- 
tures originating in propagating spinous and antisoliton 
pairs. The suppression of these features is attributed to 
the "quasi-fermionic" behavior of antisolitons with equal 
momenta. Finally we explain the different findings in 
the spectral function and the LDOS. In Sec. |V] we dis- 
cuss the effect of possible boundary bound states on the 
LDOS and in Sec. I VII we compare our results to the half- 
filled Mott insulator studie d^°'^^ previously. Finally, in 
Sec. IVIII we set our results in the context of other one- 
dimensional systems and discuss implications for STM 
experiments on quasi-one-dimensional materials. Read- 
ers mainly interested in the qualitative features of the 
LDOS may start with this section and study the detailed 
results afterwards. Technical details have been moved to 
the appendix. 



II. THE MODEL 

As the underlying microscopic model for the quarter- 
filled Mott insulator one may start with the extended 
Hubbard modeli^ 
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where ^ and cj^a are electron creation and annihilation 

operators at site j with spin a =t,i, nj,a = c] a-Cj,a, and 
Uj = Uj^-f + Uj^^. At quarter filling and for sufficiently 
large on-site and nearest-neighbor repulsions U,V > At 
this model possesses a Mott insulating phasei^ which can 
be thought of as a microscopic realization of the field- 
theoretical system we will study below. The extended 
Hubbard model (P) can be used as an effective model 
for the description of Bechgaard salts^fi; its spectral and 
transport properties have been investigated in detail us- 
ing various technique o^^'^^ . 

Instead of working with the microscopic model ^ we 
will study here the corresponding low-energy field theory. 
The continuum description is obtained by focusing on 
the degrees of freedom around the Fermi points itfep- 
We introduce slowly varying right- and left-moving Fermi 
fields as 



i?,(x)+e-*^-^"i,(x), (2) 



and similarly for the electron creation operators. Here 
flo denotes the lattice spacing, x = jao, and the Fermi 
momentum is given by fcp = 7r/4ao at quarter filling. The 
impurity potential is assumed to be strong such that we 
can model it by a boundary condition on the continuum 
electron field 



= 0) = 0. 



(3) 



Using standard bosonization of the right- and left-moving 
Fermi fields it was show n^'^'^^ that double Umklapp pro- 
cesses result in a cosine term in the charge sector, which 
is relevant at low energies and generates a gap. Thus the 
effective low-energy Hamiltonian is given by2£ 
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where, as a consequence of ([3]), the canonical Bose fields 
are subject to the hard-wall boundary conditions 

$c,s(a^ = 0) = 0. (7) 

The fields 9c, s are dual to $c,s- The charge and spin ve- 
locities Vc,s, the Luttinger parameter /?, and the coupling 
constant Qc are functions of the hopping and interactions 
in an underlying microscopic model. Within the field the- 
ory they can be viewed as phenomenological parameters. 
Experimental estimates^i for the Luttinger parameter in 
the Bechgaard salts yield fP « 0.9. 

In the bosonization procedure leading to the effective 
low-energy Hamiltonian (j4]) the explicit relation between 
the Fermi and the Bose fields is given byii>2^ 
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where the Klein factors rja- satisfy the anticommutation 
relations {ryo-,??^} = 2<5<t<t' and = 1 = -f^. 

The charge excitations are described by the sine- 
Gordon model on the half-line ([5]). For (3 < 1 the cosine 
term is relevant and generates a gap. The excitations 
are massive solitons and antisolitons which carry charges 
±e/2 respectively. They possess a relativistic dispersion 
relation; we parametrize their energy and momentum in 
the usual way using the rapidity 9 by 



E^^Acoshe*, P=— sinh6l. 
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The soliton mass A is a function^^ of the bare parame- 
ters (7c and jS but will be viewed here as a phenomenolog- 
ical parameter replacing the couphng constant g^- In the 
regime 1/2 < sohtons and antisolitons are the only ex- 
citations in the charge sector; for 0^ < 1/2 propagating 
breather (soliton-antisohton) bound states exist as well. 
At the Luther-Emery point (LEP) /S^ = 1/2 the charge 
sector is equivalent to a free massive Dirac theory^^. The 
sine-Gordon model ([5]) with the boundary condition ([7]) 
is integrabl o^^'^^ , which we will use below to calculate the 
necessary correlation functions of right- and left-moving 
Fermi fields. The spin sector © describes massless rela- 
tivistic spinon excitations which propagate with velocity 
Vs and carry spin ±1/2. We have already assumed spin 
rotational invariance, i.e. a possible Luttinger parameter 
in the spin sector is set to one. 

We note that the half-filled Mott insulator studied 
previousl y^°'^^ possesses the same effective low-energy 
Hamiltonian ([I])-®. The difference to the quarter-filled 
case studied in this article is given^^ by the bosonization 
relations ^ and (O for the right- and left-moving Fermi 
fields. This leads to differing Green functions as well as 
LDOS. We will proceed with the presentation of the re- 
sults for the quarter-filled Mott insulator and present a 
detailed comparison to the half-filled system in Sec. |Vll 



III. GREEN FUNCTION 

In order to derive the LDOS below we calculate the 
time-ordered Green function in Euclidean space, 

(r, xi, X2) = - (Obi Tr *.(r, a;i) (0, 2:2) |0b) , 

(11) 

where |0b) is the ground state of (U) in the presence of 
the boundary and t — it denotes imaginary time. At low 
energies we can linearize around the Fermi points ifep, 
which results in 
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where e.g. G^i; = - (Ob| i?.(T, xi) ^,(0, X2) |0b). For 
the calculation of the LDOS we have to set x\ = X2 be- 
low. In particular, we will focus on the spatial Fourier 
transform of the LDOS as physical properties can be 
more easily identified. The four terms in the decomposi- 
tion (I12[) then contribute in different regions in momen- 
tum space, i.e. G^^, and G^^, contribute for Q ±2fcF 
while G^^; and G^^, contribute for Q 0. In the transla- 
tionally invariant system right- and left-moving fields in 
the spin sector decouple which results in G^^, — G^^, ~ 
0. The presence of a boundary couples right and left 
sectors and concomitantly the Fourier transform of the 
Green function (jl2p acquires a non-zero component at 
Q ft! ±2fcF- We will focus on the 2fcF-part of the Green 
function and the LDOS in the following, since this pro- 
vides a particularly clean way of investigating boundary 



effects. We note that the other components of the Green 
function can be calculated in the same way. 

Due to the spin-charge separated Hamiltonian ^ the 
Green function G^^, factorizes into a product of correla- 
tion functions in the spin and charge sectors. The correla- 
tion functions in the spin sector can be straightforwardly 
derived using standard methods like boundary conformal 
field theory?!^ . On the other hand, the integrability of the 
sine-Gordon model on the half-line (jH) enables us to cal- 
culate correlation functions in the charge sector using the 
boundary state formalism^^ together with a form-factor 
expansion^Si^SiiSS. This yields the following result for the 
2fcF-component of the Green function 
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The infinite series originates from the applied form-factor 
expansion, which constitutes an expansion of correlation 
functions in the number of contributing solitons and an- 
tisolitons. In App. [A] we derive explicit expressions for 
the first three terms Gk which correspond to two-particle 
processes in the charge sector. We note that is valid 
at zero temperature. The effect of small temperatures 
T <C A on the gapless spin sector can be incorporated 
using conformal field theory^! as discussed for the half- 
filled Mott insulator in Ref. [111. In the next section we 
use the result ([T^ to derive the Fourier transform of the 
LDOS. 



IV. LOCAL DENSITY OF STATES 

The LDOS can now be calculated directly from the 
Green function. In order to analyze the physical prop- 
erties it is useful to consider the spatial Fourier trans- 
form of the LDOS as features like dispersing quasipar- 
ticles can be more easily identified. This technique was 
previously applied to the LDOS of Luttinger liquids^"— as 
weh as CDW states and half-filled Mott insulators^^. 
The Fourier transform of the LDOS for positive energies 
-E > is given by'^" 
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Here the Green function has been analytically continued 
to real times and we have taken the limit xi ^ X2 = x. 
We note that due to the assumed spin rotational invari- 
ance the LDOS does not depend on c. As mentioned 
before, we will concentrate on the 2fcF-component as it 
vanishes in the absence of the boundary and hence of- 
fers a particularly clean way of investigating boundary 
effects. For Q sa 2fcF only G^^ contributes and starting 
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from (I13|) we arrive at our main result {\q\ ^ 2kp) 



N{E, 2kF + q) = -e{E - 2A) J2 NkiE, 2kF + q), (15) 
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as well as go = q, qi = q — — sinh 0i , and (72 = 
g— 2A sinh 6*,;. Explicit integral representations for the 
scattering matrix So{0), the boundary reflection matrix 
K{0), and the function G{0) are given in App. [BJ Phys- 
ically So{0) describes the scattering of two antisolitons 
with relative momentum corresponding to the rapidity 0, 
K{0) is the reflection amplitude of an antisoliton with ra- 
pidity off the boundary [and thus incorporates the infor- 
mation on the boundary condition ([T])], and G{0) appears 
in the matrix element [i.e. form factor] of the charge part 
of ^ and with two-antisoliton states [see (|A4|]. The 
parameter ^ is defined as ^ — /3^/(l — /5')- The normal- 
ization constant-- Z2 has dimension Z2 oc j 
thus in order to obtain dimensionless quantities we mul- 
tiply the LDOS by A^/^-ZSVie-i/Z^'/y^ in ah figures. 
The singularity of the integrands in (IT51) is smeared out 
by taking 5 small but finite; we choose 5 — 0.01 unless 
stated otherwise. This mimics the broadening of singu- 
larities in experiments due to the instrumental resolution 
and the finite temperature. The LDOS vanishes for 
E < 2A as the fractionalization of an electron creates at 
least two antisolitons. 

The three terms (fTB)) constitute the leading contribu- 
tion from the form-factor expansion. As can be seen from 
the double integrals they all originate in two-particle pro- 
cesses in the charge sector. The number of boundary 
reflection matrices K{0) indicates that Nk involves k in- 
teractions with the boundary, i.e. the terms (|16p contain 
no, one, or two reflections of charge excitations. The sub- 
leading contributions involve at least three particles in 
the charge sector which either come from a higher num- 
ber of particles in the intermediate state or from higher- 
order processes due to the boundary. Similar terms were 




FIG. 1: Constant energy scan for E = 3A: \N{E, 2kF+g)\ and 
ArgA''(£;,2fcF +q) for Vc = 1.2 and = 0.9 (upper panel) 
and = 0.7 (lower panel). Both plots are shown on the same 
scale. W e observe a peak at g = and dispersing features at 
q = ±2^{E - A)2 - as well as q = ±2^/ E'^ - AA?/vc. 

For q <Q the dispersing features are strongly suppressed. 



thoroughly analyzed for the LDOS in a CDW state at the 
LEPl''^ as well as the spectral function in the Ising model 
with a boundary magnetic field"^^ j they were found to be 
negligible in both cases. We further note that the sup- 
pression of sub-leading terms in the form-factor expan- 
sion for bulk two-point functions is a well-known feature 
of massive integrable field theories^!^. We thus expect 
the higher-order corrections to ([15]) to be negligible in 
the low-energy regime we consider here. 

The constant energy scan of the Fourier transform of 
the LDOS (flSl) is shown in Fig. [1] We observe a singu- 
larity at g = 0, i.e. Q = 2A:f, which is indicative of the 
pmnmg of the 2fcF-SDW at the boundary (see Sec. lIV A|l . 
Furthermore there is a peak at g > 0. The constant 
momentum scans shown in Fig. [2] reveal that this peak 
shows dispersing behavior and splits above a critical mo- 
mentum q > qo [qo will be defined in Eq. (j22p ]. More pre- 
cisely these features follow the dispersion relations p9)) 
and ((20|) respectively (see Sec. I IV Bp . The existence of 
dispersing features is at first surprising since the spectral 
functiortii of the translationally invariant system does 
not show any dispersing quasiparticles. The processes 
contributing to the LDOS (fT5|) can be thought of as aris- 
ing from creating an electron at a position x, the fraction- 
alization of this electron into at least one spinon and two 
antisolitons, the propagation of these elementary excita- 
tions through the system, and their subsequent recombi- 
nation at the position x. This allows a natural interpre- 
tation of the leading terms in the form-factor expansion 
(|16p in terms of propagating antisolitons and spinous. In 
the next sub-sections we will analyze the different terms 
(fT6|) separately. Afterwards we comment on the different 
findings in the spectral function and the LDOS. 
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FIG. 2; \N{E,2kF + q)\ for = 0.9 and Vc = 1.2 w^. The 
curves are constant g-scans which have been offset along the 
y-axis by a constant with respect to one another. The LDOS 
is dominated by the peak at g = 0. We further observe dis- 
persing features at Eic = A + + {vcq/2)^ as well as 
Eica = A + Ay^l - {vs/v^Y + for g > go. 

A. No reflections in the charge sector: 
NoiE,2kF +q) 

Let us start with the analysis of the first term in the 
form- factor expansion ([T5| , which is shown in Fig. [3l We 
first note that this term is dominated by a singularity 
at g = 0, i.e. Q = 2fcF- This singularity is due to the 
pinning of the 2fcF-SDW at the boundary. As fcp — 7r/4ao 
we find a spatial modulation with periodicity 27r/2fcp = 
4ao as sketchedi^i^ in Fig. |31b. Close to g = the first 
term in the form-factor expansion and thus the whole 
LDOS behaves as^ 

N(E,2kF + q) ^ No{E,2kF + q) (17) 

We note that the exponent is independent of the Lut- 
tinger parameter /3 (see Fig. |3la). This singularity is 
similar to the ones observed in the LDOS of Luttinger 
liquida^i^ as well as CDW state O' S' ""' where, however, the 
exponents depend on the interaction strength and thus 
the Luttinger parameter. 

Furthermore, Nq possesses a weak dispersing feature 
following (see^^ Fig. Oc) 

EM = 2A + ^, q>0. (18) 

As already discussed, the created electron fractionalizes 
into two antisolitons and (at least) one spinon which can 
propagate through the system. In this way of thinking 
the dispersing feature at Eg arises from two antisolitons 
with zero momentum (solid lines in Fig.[3ld) contributing 
an energy 2 A and a massless soliton with momentum q 
(dashed line in Fig. Od). The appearance of Vs/2 in Eg 
is due to the fact that the spinon has to propagate to 
the boundary and back, thus covering the distance 2x in 
time t. 



a) 




FIG. 3: (Color online) No{E,2kF + g) for «c = l.6v,. a) 
Logarithmic plot close to the singularity at g = for S = 
0.001. The dotted line is a guide to the eye; it has gradient 
-1/2. b) Sketcbi^i^ of the 2fcF-SDW pinned at the boundary, 
c) \No{E,2kF + q)\ for = 0.9 and different momenta. We 
observe a dispersing feature at Es = 2A -I- Vsq/2 indicated 
by small arrows, d) Sketch of the process giving rise to the 
propagating feature. The electron decomposes into two zero- 
momentum antisolitons (solid lines) and one spinon (dashed 
line) reflected at the boundary. 



B. One reflection in the charge sector: 
Ni(E,2kF + q) 

The next term in the form-factor expansion (jlSp con- 
tains one boundary reflection matrix K(9). Thus it is 
natural to interpret it as arising from processes involv- 
ing the reflection of one antisoliton at the boundary. In 
Fig. IHa we show Ni for different momenta. We clearly 
observe a dispersing feature (indicated by up-arrows) fol- 
lowing 

i?ic(g) = A + yA2 + (M)'. (19) 

This feature originates in the process where one anti- 
soliton propagates with momentum q while the other 
antisoliton and the spinon possess zero momentum (see 
Fig. Hlb). The first term in Eic is the rest mass of the 
zero-momentum antisoliton while the second term is the 
energy of the propagating antisoliton. 

Furthermore, when q exceeds the critical value go a 
second dispersing feature appears (indicated by down- 
arrows in Fig. Hla) at 



i?ies(9) = A-KAyi-(^JJi)%M. (20) 

The critical momentum q^ can be determined from the 
condition that the spin and charge excitations have the 
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FIG. 4: (Color online) a) \Ni{E,2kF + q)\ for fi^ = 0.9, = 
1.6 Va, and different momenta. We observe dispersing features 
at Eic = A + A'^ + {vcq/2Y indicated by up-arrows and 
Bics = A+Aiyi — {ys/vc)'^-\-Vsq/2 indicated by down-arrows, 
b) Process resulting in Eic- The electron decomposes into one 
antisoliton and spinon with zero momentum as well as one 
antisoliton reflected at the boundary, c) Process resulting in 

-Blca. 



same group velocity 



which leads to 



dq 



2 



2Ai 



(21) 



(22) 



Thus this feature can be thought of as arising from an 
antisoliton with momentum qo and a spinon carrying mo- 
mentum (7 — go, while the second antisoliton stays at rest 
(see Fig. Hie). The fact that the propagating antisoliton 
and the spinon have the same group velocity results in 
an increased probability that both simultaneously arrive 
at position x after reflection at the boundary. 

The observed splitting of the quasiparticle peak is rem- 
iniscent of what is found for the single-particle spectral 
function in the translationally invariant, half- filled Mott 
insulatorj^ as well as the LDOS of a CDW state in the 
presence of a boundary in the regime of attractive inter- 
actionsiSiii. The peak splitting is a consequence of the 
curvature of the antisoliton dispersion (and thus of the 
charge gap) together with the condition Vc > Vg- Hence it 
cannot be observed in the Luttinger liquid case^"— where 
both sectors are massless. 



C. Two reflections in the charge sector: 
N2iE,2kF +q) 

The leading term in second order in the boundary re- 
flection matrix is A^2- As in the previous term we find 



FIG. 5: (Color online) a) \N2{E,2kF + q)\ for = 0.9, 
Vc — 1.6 Us, and different momenta. We observe dispers- 
ing features at E2c ~ 2^ + (fcfJ'/4)^ (up-arrows) and 
E2ca = 2A Vl - (vs/vcy + V sq/2 (down-arrows). Processes 
resulting in (b) E2c and in (c) £'2cs. 



two dispersing features. First, indicated by up-arrows in 
Fig. Oa, there is a feature following 



E2ciq) 



(23) 



This feature can be thought of as arising from the prop- 
agation of both antisolitons with momenta q/2 and thus 
the same velocity, while the spinon has zero momentum 
(see Fig. [S]b). The second feature is found at 



E,.4q)^2A\ll-['-^) +M, 



(24) 



provided q exceeds the critical momentum 2qo. This fea- 
ture is indicated by down-arrows in Fig. Oa. The in- 
terpretation is that both antisolitons carry momentum 
qo while the excess momentum q — 2qo is carried by the 
spinon (see Fig.Elc). 

We note that at finite momentum q > the term A^i 
dominates the LDOS as can be seen by comparing the 
scales in Figs. [He, Hja, and [Ha, respectively. More im- 
portantly, the dispersing features originating in A'^i are 
much more pronounced than those resulting from Nq and 
N2. We attribute this to the fact that, at least in the 
naive interpretation of Figs. [Db and c, the antisolitons 
involved in processes contributing to A^i possess different 
momenta as one of them is at rest while the other propa- 
gates through the system. In contrast, the dispersing fea- 
tures in iVo and N2 originate from processes in which the 
two antisolitons possess equal momenta (see Figs. [3|d as 
well as[5lb and c, respectively). Thus their relative rapid- 
ity vanishes and the scattering matrix is So{9 = 0) — —1 
(see App. [BJ. The resulting "quasi- fermionic" behavior 
leads to a strong suppression of the dispersing features 
in iVo and N2 via the Pauli principle. 
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Finally let us comment on the different findings in the 
spectral function--^ of the translationally invariant sys- 
tem and the LDOS (fT5|l in the presence of a bound- 
ary. The spectral function is obtained from the Green 
function G{t,x,0) via Fourier transformation with re- 
spect to t = —IT and X. The contributing processes 
require the propagation of all excitations, and in par- 
ticular both antisolitons, from (0,0) to {t,x). Thus one 
may expect a dispersing feature at i? = 2-\/A2 + {vcq/2y 
from the antisolitons propagating with the same mo- 
mentum q/2. In addition, above the critical momentum 
2(7o — '2,/^Vs/vc\/v^—vi a linearly dispersing feature fol- 
lowing E = 2A^1 — (us/wc)^ + Vsq is expected. Indeed 
these are exactly the conditions™ for the threshold be- 
low which the spectral function vanishes. Above the 
threshold the spectral function is rather featureless, in 
particular there are no dispersing peaks associated with 
the propagation of charge or spin degrees of freedom. 
The absence of dispersing features may be attributed to 
the "quasi-fermionic" behavior of antisolitons, as all pro- 
cesses contributing to the spectral function will contain 
antisolitons with equal momenta. In contrast, the pro- 
cesses contributing to the LDOS dU]) are given by the 
fractionalization of an electron at position x, the propa- 
gation of the two antisolitons and the spinon through the 
system, and the subsequent recombination at the original 
position X. In particular, processes where the momentum 
q is solely carried by one of the antisolitons or a spinon- 
antisoliton pair will contribute, resulting in the dispersing 
features at Eic and i?ics, respectively. The observation 
of these features thus underlines that the spectral func- 
tion and the LDOS probe rather different physical pro- 
cesses. While the spectral function requires the propaga- 
tion of all excitations between two points in the system, 
the LDOS also probes processes where only one or two 
excitations propagate to the impurity and back. 



V. LOCAL DENSITY OF STATES IN THE 
PRESENCE OF BOUNDARY BOUND STATES 

Up to now we have concentrated on the boundary con- 
ditions ([3|) for the electron field or equivalently (O for 
the canonical Bose fields. In this section we will consider 
more general phase shifts for the electron fields which 
lead to arbitrary Dirichlet boundary conditions for the 
Bose fields. As is well known, such boundary condi- 
tions can result in the existence of boundary bound states 
(BBS's), which are expectedi^iii to manifest themselves 
in the LDOS as features within the spectral gap. Here 
we will focus on general Dirichlet boundary conditions in 
the charge sector 



$c(x = 0) = 



13- ^-/?' 



(25) 



while in the spin sector we still assume $3(2; = 0) = 0. 
In the boundary state formalism^^ applied here the infor- 
mation about the boundary condition (1251) is encoded in 




FIG. 6: \N{E,2kY+q)\ for = 0.5, $° = -4, imAvc = l.2vs. 
The curves are constant g-scans which have been offset along 
the y-axis by a constant with respect to one another. The 
existence of a BBS results in spectral weight within the two- 
soliton gap E < 2A. 



the boundary reflection matrix K(6). In particular, the 
existence of a BBS will show up as a pole in the physical 
strip < lm9 < 7r/2. As was shown in Ref. [H this will 
be the case provided 



(26) 



the pole at = i7 corresponds to a BBS with energy 



>bbs 



A sin 7, 7 



2 - 2/?2 • 



(27) 



When calculating the Green function the pole of the 
boundary reflection matrix will lead to two additional 
terms, and Q4, in the form- factor expansion (|13p . This 
in turn yields two additional terms in the LDOS (I15p . 

To be explicit, let us consider the LEP (3"^ — 1/2 in 
the following. The three leading terms in the form-factor 
expansion of the LDOS are still given by (|16|) where the 
boundary reflection matrix now depends explicitly on 



K{e) 



sin (if 



V8' 



^2 + vl) 



(28) 



This obviously possesses a pole in the physical strip < 
lm9 < 7r/2 provided $[! < — 7r/\/2 corresponding to the 
existence of a BBS. The resulting two additional terms 
in the LDOS read (k = 3, 4) 

^0 g-i(77-V2<l?^)/8 



Nk{E, 2fcF +q) = Z2y/U: cos 

de e{E - Acoshe - Ehhs) 

2^ V-E- Acosh6'-£;bbs 
hk{0) 



r3/2 



Vsqk - 2(E - A cosh 6* - E\ 



bbs 



IWsKbbs 



(29) 
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where 



sinh 



+ Utt + V2<S>°) 



hi{e) = -K{e + i^)e'^'^ sinh^ 



bbs 



-A sin 



$0 



'^bbs 



^^-i(7r + ^/2<i>?) 



2A $" 
cos—, 



as well as 93 = g and <?4 = (? — ^ sinhfl. We note that 
Kbbs > has the dimension of an inverse length. The 
spatial dependence of the terms Qz and (see App. \Q 



is dominated by 



(x < 0), thus it is natural to 



interpret l/^bbs as the width of the BBS. 

The terms ([^ possess a threshold at i? = A + E^bbs 
corresponding to the creation of the boundary bound 
state and an additional antisoliton. Interestingly, al- 
though A^3 and N/^ separately possess jumps at the 
threshold, the sum of both terms is continuous. This is 
reminiscent of the BBS contributions in the Ising model 
with a boundary magnetic field^^, where a cancellation 
of singularities yields a smooth spectral function at the 
threshold. 

The LDOS in the presence of a BBS is shown in Fig. [51 
We clearly observe a contribution within the gap E < 
2A. In contrast to BBS's in the CDW statei^iii there is, 
however, no singularity at E = E'bbs as the underlying 
processes involve the creation of a spinon, an antisoliton, 
and the BBS. As before the LDOS is dominated by a 
singularity at g = and we observe dispersing features 
following ([T9| and pO)) respectively. 



VI. 



COMPARISON TO HALF-FILLED MOTT 
INSULATOR 



The LDOS of a half-filled Mott insulator in the pres- 
ence of a boundary was analyzed in detail^ in Ref. [111. 
As already mentioned, the effective low-energy Hamilto- 
nian is also given by (jl])-® but the expressions for the 
right- and left-moving Fermi fields in terms of the Bose 
fields differ from ([5]) and From the physical point 
of view right- and left-moving Fermi fields create and 
annihilate individual antisolitons in the half-filled case 
whereas they create and annihilate antisoliton pairs in 
the quarter-filled system. Thus the gap in the LDOS is 
given by A in the former and 2A in the latter system. 
Apart from this we find the following similarities and dif- 
ferences. 

In both cases the LDOS possesses a singularity at 
q = which behaves as N{E, 2kF + q) ~ l/^/U^ indepen- 
dently^ of the Luttinger parameter j3. This singularity 
is due to the pinning of the 2fcF-SDW at the boundary, 
which possesses a wave length 2ao or 4ao respectively. 

Beside this singularity the LDOS of the half-filled 
M ott insulator shows dispersing features at E = 
a/A^ -I- (wcg/2)2, at E ^ A + Vsq/2, and, for q > qo, 
at E = Ay/1 - iy^lv^f + q/2. An interpretation is 



obtained by noting that the electron decomposes into 
one spinon and one antisoliton, which then propagate 
through the system. As the leading processes involve only 
individual antisolitons, the " quasi-fermionic" behavior of 
the antisolitons discussed at the end of Sec. lIV Cl does not 
affect the LDOS. In contrast, in the quarter-filled system 
we observe only two dispersing modes following (1191) and 
([20]) respectively (see Fig. [5]) . All other dispersing fea- 
tures discussed above are strongly suppressed due to the 
"quasi-fermionic" behavior of the antisolitons. 

Finally, in the presence of a BBS with energy E'bbs < 
A the LDOS of the half-filled Mott insulator pos- 
sesses a non-dispersing singularity with N{E, 2fcF + q) ^ 
tj \JE — Sbbs which originates from the creation of a 
spinon and the BBS. In contrast, in the quarter-filled 
case the LDOS is smooth at the threshold = A -I- i?bbs 
which originates from the simultaneous creation of one 
antisoliton and the BBS. 



VII. DISCUSSION AND EXPERIMENTAL 
SIGNATURES 

In this section we compare qualitative aspects of the 
LDOS of various one-dimensional systems in the presence 
of a boundary or strong impurity potential. We discuss 
the spatial Fourier transform N{E,Q) defined in (|14p . 
Specifically we consider Luttinger liquids^, half-filled 
Mott insulators (or CDW states and the quarter- 
filled Mott insulator investigated in Sec. IIVI We focus on 
the Q ~ 2fcF-component of the LDOS as it vanishes in 
the translationally invariant systems and thus provides 
a particularly clean way to investigate boundary effects. 
(We note that the LDOS in the small-momentum regime 
Q « behaves qualitatively similar.) The first feature, 
which all three systems have in common, is a singularity 
of the LDOS N{E, 2kp + q) at q — 0, which originates 
in the pinning of a charge or spin density wave at the 
impurity. In the case of a Luttinger liquid the strength 
of this singularity depends on the Luttinger parameter 
Kc and thus on the interactions between the electrons, 
whereas in both the half-filled and quarter-filled Mott in- 
sulator the LDOS behaves as^ N{E,2kF + g) ~ 1/^- 
Furthermore, as was pointed out in Ref. IH, the gap of 
the LDOS in the half-filled Mott insulator is given by the 
soliton mass A, which equals the thermal activation gap 
and is half of the gap observed in optical measurements. 
In contrast, in the quarter-filled system the LDOS pos- 
sesses a gap of 2 A, which is equal to the optical gap but 
twice the thermal gap. 

In addition, the LDOS possesses various dispersing 
features (see Fig. [7|). In the Luttinger liquid case one 
observes two linearly dispersing modes corresponding to 
individually propagating spin and charge degrees of free- 
dom respectively (see Fig. [71a). In the half-filled Mott in- 
sulator an electron decomposes into massless spin excita- 
tions and one charge excitation with mass A. Propagat- 
ing spin excitations give rise to a linearly dispersing mode 
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Luttinger liquid half-filled MI quarter-filled MI 

FIG. 7: Schematic picture of the dispersing features observ- 
able in the spatial Fourier transform N{E, 2fcp + q) of the 
LDOS in different one-dimensional systems, a) In a Luttinger 
liquid one finds two linearly dispersing modes, b) In a half- 
filled Mott insulator (MI) there exists a gap A. One observes 
one linear and one curved dispersion. From the latter a third 
feature splits off at 5 = go- c) The quarter-filled Mott insu- 
lator possesses a gap 2A. One further observes one curved 
dispersion from which a linear mode splits off at q = qq. In 
both (b) and (c) the thermal activation gap is given by A. 



with velocity Wg, while the propagating charge excitation 
results in a curved dispersion relation. In addition, there 
exists a critical momentum q = qq ed, which the group ve- 
locity of the charge excitation equals Vg- At this momen- 
tum a third linearly dispersing feature splits away from 
the charge mode, which has its origin in the collective 
propagation of spin and charge degrees of freedom with 
equal velocities (see Fig. Hb). We note that the qualita- 
tive behavior discussed for the half-filled Mott insulator 
is also expected in other systems possessing gapless exci- 
tations in one sector and gapped ones in the other like, for 
instance, CDW states. Finally, in the quarter-filled Mott 
insulator an electron decomposes into massless spin exci- 
tations and two charge excitations with masses A. This 
results in a curved dispersion originating in the propaga- 
tion of one charge excitation while the second one stays 
at the position of the STM tip. This dispersion splits at 
the critical momentum qq due to the collective propaga- 
tion of one charge excitation and additional spin degrees 
of freedom (see Fig. He). We note that one does not 
observe individually propagating spin excitations. The 
relevant processes would require the two charge excita- 
tions to possess equal momenta, which is forbidden by 
their " quasi- fermionic" behavior (see Sec. IIV Ct . 

To sum up, at sufficiently large momenta q we observe 
two linearly dispersing modes in a Luttinger liquid, two 
linear and one curved dispersion relation in a half-filled 
Mott insulator, and one linear and one curved dispersion 
in a quarter-filled Mott insulator. This clearly demon- 
strates that characteristic properties of the bulk state of 
matter and its electronic excitations can be infered from 
studying spatial modulations of the LDOS in the pres- 
ence of an impurity. 

The LDOS discussed in this article is directly related 
to the tunneling current measured in STM experiments. 
If the density of states in the STM tip is approximately 
independent of energy, the tunneling conductance will be 
given by the thermally smeared LDOS of the sample at 



the position of the tipi 

^l^^^ J dEr{E~eV)N{E,x), (30) 

where f{E) denotes the Fermi function. N{E,Q) can 
then easily be obtained via spatial Fourier transforma- 
tion. In this context the models discussed above ap- 
ply to quasi-one-dimensional materials at energies above 
the cross-over scale to three-dimensional behavior. This 
situation might be experimentally realized for example 
in stripe phases of high-temperature superconductors^'^, 
carbon nanotubes'^-"*^, Bechgaard and Fabre salts^^, self- 
organized atomic gold chains'*^, and two- leg ladder ma- 
terials^-. 

As discussed above already the knowledge of qualita- 
tive properties such as the presence of a spectral gap and 
its relation to the thermal activation gap as well as the 
number and positions of singularities can reveal substan- 
tial information about the electronic properties and ex- 
citations of quasi-one-dimensional materials. A more de- 
tailed analysis of quantitative features like the strength of 
singularities or the precise positions of dispersing modes 
can even lead to estimates for the Luttinger parameter 
and the velocities of charge and spin excitations. There 
are, of course, experimental limitations. Obviously the 
experimental resolution and thermal broadenings have 
to be sufficiently small to clearly distinguish the differ- 
ent singularities and dispersing modes. Moreover, the 
simplyfied relation (150)) does not contain^ the matrix el- 
ements for the tunneling from the sample into the tip, 
which may introduce a non-trivial spatial dependence of 
the tunneling conductance on the LDOS and thus ham- 
per a natural interpretation in terms of propagating ex- 
citations. 



VIII. CONCLUSIONS 

In conclusion, we have studied a quarter-filled Mott in- 
sulator in the presence of an impurity, which was modeled 
by a boundary condition for the electron fields. In this 
system we calculated the 2A:F-component of the spatial 
Fourier transform of the LDOS using the boundary state 
formalism together with a form-factor expansion. The 
LDOS is dominated by a singularity at Q = 2k-p, which 
is indicative of a pinning of a SDW at the impurity. Fur- 
thermore, we observed several dispersing features above 
the two-soliton gap E < 2A. We identified their physi- 
cal origin as spin and charge degrees of freedom, which 
propagate between the STM tip and the impurity, and at- 
tributed their relative strength to the "quasi-fermionic" 
behavior of charge excitations with equal momenta. This 
also explains the absence of dispersing features in the 
spectral functionF^ of the translationally invariant sys- 
tem and underlines that the spectral function and the 
LDOS probe rather different physical processes. We fur- 
ther showed that a BBS in the charge sector leads to 
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a non-vanishing LDOS within the two-sohton gap. Fi- 
nahy, we discussed our results in the context of other 
one-dimensional systems, i.e. Luttinger liquids and half- 
filled Mott insulators, and argued that the measurement 
of the spatial modulations in the LDOS can be used to 
extract detailed informations about the electronic states 
of quasi-one-dimensional systems. 
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Appendix A: Calculation of the Green function 

In this appendix we derive the leading terms Qk in the 
form-factor expansion of the charge part of the Green 
hmction (|13p . i.e. we determine the correlation function 

(Ofl(r,a:i)Ol(0,X2)) (Al) 



in the sine-Gordon model ([5]) on the half line. Here the 

operators 0|j and are given by the charge parts of 
(HI) and Q respectively. We calculate (|Aip using the 
boundary state formalism-^ together with a form-factor 
expansionS^iSSi^. The same technique was previously ap- 
pliediSiiii to determine the Green function of a CDW 
state in the presence of an impurity. Here we restrict 
ourselves to the main steps of the derivation and con- 
centrate on the differences as compared to Ref. [ll|, to 
which we refer implicitly for all details and definitions 
not discussed here. 

By performing a rotation in Euclidean space we map^ 
the problem onto the sine-Gordon model on the real axis. 
The boundary condition is translated into an initial con- 
dition encoded in a boundary state. Expanding the latter 
in powers of the boundary reflection matrix and inserting 
a resolution of the identity between the operators in (|A1[) 
we find 

oo oo 

2m(T, a;i,X2), (A2) 

n—O m—0 

where we have defined the auxiliary functions 



Cn2miT,Xi,X2) 



'2™ to! n! (277)™ 



dei...de„ ^ai6i(^,)_ _^a„fc,„(g^) 



X {0\OR{T,xi)K,...,di)^ 



^{e^,...,e,,\o{[Q,x2)\~e[,e[, 



(A3) 



■•I ^rm^ml ai,bi,...,am.b„ 



Here the states |6'„, . . . , (and so on) are scatter- 

ing states of solitons and antisolitons, the indices ai,bi, Ci 
take the values =f1 and label the corresponding U(l) 
charges, and energy and momentum of solitons and an- 
tisolitons are given in terms of the rapidities 9i and 9'^ 
by (fTU]) . We note that (jA2p constitutes an expansion in 
two parameters: (i) the number of particles (with rapidi- 
ties 6i) in the intermediate state and (ii) the number of 
reflections of particles at the boundary [the number of 
boundary reflection matrices K{9'j)]. 

The scattering of solitons or antisolitons with rapidi- 
ties 9i and 6*2 is encoded in the scattering matrijs^ 
l^i>^2),, = ^^^(^1 - 92) \92,9i)^^. The boundary reflec- 
tion matrixes. K'^^{9) similarly describes the reflection at 
the boundary, \9)^ = ii'"''(if - 9) |-6')j, where a = ± for 
a = =F. For Dirichlet boundary conditions one has^i^ 
K++{9) = K—{9) = 0. 

Furthermore, 0{t,x) = e~^^ e~^^'^Oe^^'^e'^^ , where 
H and P are the Hamiltonian and the total momentum 
of the sine-Gordon model on the infinite line. The oper- 
ator Or {O]^) creates (annihilates) two antisolitons, i.e. 
increases (decreases) the U(l) charge by 2. The necessary 



form factors were derived by Lukyanov and Zamolod- 
chikov'^^. In our conventions they are given by 

{0\Or^l\9,,92)__ 

- V^e±'-/«e±(«^+«^)/«G(0i-02), 

where the normalization constant Z2 and the auxiliary 
function G{9) are stated in App. [B] The form factors 
satisfy the form-factor axioms as stated in Ref. [ll| with 
Lorentz spin s{Or/l) = ±1/4 and semi-locality factor 
1_{Or/l) — e^"^/*. In the following we will evaluate the 
terms C20, C22, and C24 which result in Qq, Qi, and G2, 
respectively. 

1. Derivation of Qo 

The first term to be evaluated is C20 ■ It does not con- 
tain any information about the boundary in the charge 
sector, i.e. it contains no boundary reflection matrices. 
A similar term was calculated in Ref. [l3 in the deriva- 
tion of the spectral function in the translationally invari- 
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ant system. After shifting the contour of integration, 
&i,2 &1.2 + i7i'/2, we directly obtain 



^^0 = C2 



deid02 



\Gi9i 



01 



where the center-of-mass coordinates are defined by a: = 
{xi + X2)I1 < and r = xi ~ X2 < 0- Analytical contin- 
uation T — >■ it to real times, taking the limit r — >■ 0, and 
Fourier transformation (jl4p with respect to t and x yield 
the first term No{E, 2kF + q) in the form-factor expansion 
of the LDOS (fl^l). 



Derivation of Qi 



We start with C22 defined in (|A3|) . As the operator 
Or creates two antisolitons we deduce ci = C2 = — ; 
K++{e') = K—{e') = yields b = d. The matrix el- 
ement of 0\ contains incoming and outgoing particles 
and thus possesses kinematical poles. We deal with these 
terms following Smirnov^^ and analytically continue the 
form factor as (see Ref. [Il| for a detailed discussion of 
this procedure) 



-- -(0i+iO,(?2+iO|Oi|-0',0\a 



(A5) 



+2lr5{e^ + 9') 5^, st{-2e') -(^^l \e')^ 

+27rd{02 + 0')SbeSil{^20') 

xsiU0i-e2)''{0i\o[\e')^. 

Here we have already omitted the positive imaginary 
parts of the rapidities in the form factors in the 2nd to 
5th lines as they do not possess kinematical poles. In- 
serting (jASP into C22 yields to different contributions: 
The first line leads to a term containing three integra- 
tions over rapidities and is thus a sub-leading correc- 
tion; we will not evaluate it here. The 2nd to 5th 
lines yield using the boundary cross-unitarity condi- 

tion2& K''''{e) = S'^l{20) K'^''{-e) as weU as unitarity^ 
'6162 ( 



si\l\{e)slili{-0) 



5^^ s" 



de' d0 

27r 27r 



K+-{0') {O\OR\e,0')_ 



01 



„ „iAT(sinhe+sinhe') ^ (r cosh 9+22; cosh 9') 



(A6) 



In order to proceed let us first assume — 0. This im- 
plies that the boundary reflection matrix does not depend 
on the U(l) charge, i.e. K+-{0') = K-+{0') = K{e'), 
where K{0) is explicitly stated in App. [BJ In particular, 
K{0') is analytic in the physical strip < Im6'' < tt/2. 



Now we shift the contours of integration 0^ 0' 0,0' + 
m/2 (we note that the contribution originating in the first 
term of (|A5|1 does not possess any poles in the physical 
strip < lm0i < tt/2) and apply the crossing relation in 
the second form factor, 



-(0 + if I 0[ \~0' - if) =+{-0' + if I Ol 



if): 



-i7r/4 



(O|0z 



; 37r 



if): 



where the phase factor e"''^/'* is due to the semi- 
localit y^^i^^'^^ of the operator Ol with respect to the 
fundamental fields creating the excitations. Finally we 
use (IA4I1 to obtain 



i7r/4 



d0' d0 

27r 27r 



K{0' + i§)G{0-0')G{0 + 0'y 



X g^74 gi:^(''sinhe-|-2a;sinhe') Ar (cosh e+cosh 9') 



CJi = Z2 e 



Fourier transformation p4|) directly yields the second 
term Ni{E, 2kF + q) of the LDOS. 

On the other hand, if $[! < — /^tt, the boundary re- 
fiection matrix {0') possesses a pole at 0' — i'y. At 
the LEP the residue is given by 2icos($°/y2) (see (jM])). 
Thus when shifting the contours of integration in (|A6p 
we obtain the additional BBS contribution 



2i Z2 cos 



, ^ -i(7r-\/2*«)/8 -BbbsT Kbbsa; 



d0 

2^ 



sinh 



V2 

0+^{-K + ^/2<^l) 



• sinh 9 g — Ar cosh 9 



Fourier transformation yields N^{E, 2k-p + q). 

3. Derivation of Q-^ 

Finally we calculate the leading term involving the re- 
flection of two antisolitons, which originates from C24. 
The analytic continuation of the form factor of 0\ reads 

{0l,02\O\ \~O[,0[,~'0'2,0'2)^g^hb 
= --(01+iO, 02 + iO\ 0[ \-0[,0[, -0'2, 0'2)a^bi 

+SZ+i0[ + 0'2)6+J+b 

X--{0^,02\0[,0'2}— {O\O[\-0[,-0'2) 

+S-+{-20[)SZ+{-(^[+0'2)S+b 

X — {0,,02\--0[,0'2)— (O|Ol|0'i,-C 



(A7) 



+S^;{~20'^)SZt{0[~0^2)S+a 



9' - 



\O[\-0[,0'2] 



++ 



++ 



+S-+{~20[) S7j+{-20'2) SZ+{-0[ - 0'2 



'bb 

X — {0^,02\-0[,-0'2)— {O\O[\0[,0'2)++, 

where the scalar products are evaluated using the 
Faddeev-Zamolodchikov algebra^^i^ 

-{01,021 O[,0'2)- = (27r)2 [<5(0i - 0'2) 6(02 - 0[) 

+ SZZ{0[ - 02) 6(01 - 0[) 6(02 ~0'2)]. 
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Therefore the 2nd to 5th Hne in (IA7|) yield (we also use 
the scattering axiom for the first form factor) 



1 



2.7-0O (27r)2 

x{O\On\el0',)__ {o\ol\-e[,-9',)_ (^8) 



Let us again first assume $2 = 0. If we shift the contours 
of integration 9[ — >■ 0^ + i7r/2 and use the crossing relation 
as well as a Lorentz transformation, 

(0|ol|-^^l-if,-^?^-if), , 



\OL\-e: 



21 " 



we obtain 

G2 



^ „i^/4 



X szXie'i + 6*2 + 17^) iGlf'i - 6*2)1 

X e(el+e2)/4 g2i^:c sinh 9[ ^-At cosh t 



matrix between solitons and antisolitons possesses the 
integral representation^ 



exp 



dt . 6t sinh(|(l - I))' 

— sin — ^ — 

t 7rf sinh | cosh ^ 



where ^ = — /3^). A similar integral representation 

for the boundary reflection matrix reads^ (for Dirichlet 
boundary conditions with $2 = 0) 



K(e) = -cos ( _ + i- ) i?o(i| _ e) a(if - 



i?o {0) — exp 
ct(6') =exp 



2i 



dt . et . 3t sinh(|(l- i)) 

— sm — smh ; 

t TT^ 4^ sinh J sinh j 



dt , et sinh(|(l + if)) 
2i / — sm — — r-r r-i — 

t smh t cosh |- 



The auxiliary function G{6) appearing in the form factors 
21 is given hy^ 



N2{E, 2kF +q) is now readily obtained via Fourier trans- 
formation. 

Starting from (jASp the second term Nn{E,2k-p + q) 
related to the BBS is obtained by picking up the poles of 
K+~[6[) and K+^{92) at 2 = 17 respectively. At the 
LEP the result is 



^-i(;r-\/2*°)/8 g--EbbsT gftbbsa; 



Qa = 2i Z2 cos — i e 
V2 



°"^A-(.' 
27r 



if) sinh^ 



i(7r + V2<i>?) 



Appendix B: Explicit expressions for form factors 
and scattering matrices 

In this appendix we collect some formulas on the sine- 
Gordon theory used in the present article. The scattering 



G{e) = -iCisinh- 



X exp 



exp 



dt sinh^(t(l + i|)) sinh(t(e-l)) 
t sinh(2i) cosht sinh(t^) 



dt sinh^ I sinh(t(C- 1)) 



t sinh(2i) cosht sinh(t^) 



while for the normalization factor Z2 one finds numeri- 
cally from the integral repesentation derived in Ref. |3i| 

Z2 ~ 7.7320 A1681/1440 for p2 ^ 
Z2 ~ 20.1857 Ai^4^/"20 for ^2 ^ q j 

The expressions at the LEP are obtained by setting ^ = 1 
{^^ = 1/2), i.e. So{9) = -1, K{9) = itanhf, G{e) = 
-isinhf, and Z2 w 38.5519 
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